per 


WORLD INTELLECTUAL PROPERTY ORGANIZATION 
International Bureau 



INTERNATIONAL APPLICATION PUBLISHED UNDER THE PATENT COOPERATION TREATY (PCI) 


(51) International Patent Classification 5 
H04S 1/00, H04R 29/00 
H03G5/16 


Al 


(11) International Publication Number: WO 90/00851 

(43) International Publication Date: 25 January 1990 (25.01.90) 


(21) International Application Number: PCT/GB89/00773 

(22) International Filing Date : 7 July 1989 (07.07.89) 


(30) Priority data: 
8816364.7 
8910847.6 


8 July 1988 (08.07.88) GB 
11 May 1989(11.05.89) GB 


(71) Applicant (for all designated States except US): ADAPTIVE 

CONTROL LIMITED [GB/GB]; Norwich Union 
House, 102-108 Above Bar Street, Southampton S09 
5HJ (GB). 

(72) Inventors; and 

(75) Inventors/Applicants (for US only) : NELSON, Philip, Ar- 
thur [GB/GB]; 19 Chetwynd Drive, Southampton, 
Hampshire S02 3HY (GB). ELLIOTT, Stephen, John 
[GB/GB]; 26 Elm Road, Winchester, Hampshire S022 
5 AG (GB). STOTHERS, Ian, McGregor [GB/GB]; 89 
Newtown Road, Eastleigh, Hampshire S05 4BX (GB). 


(74) Agent: BOUTLAND, J., A.; Barker, Brettell & Duncan, 
Prudential Buildings, Room 24, 97-101 Above Bar 
Street, Southampton S09 4GT (GB). 

(81) Designated States: AT (European patent), AU, BE (Euro- 
pean patent), CH (European patent), DE (European pa- 
tent), DK, FR (European patent), GB (European pa- 
tent), IT (European patent), JP, LU (European patent), 
NL (European patent), SE (European patent), US. 


Published 

With international search report 


(54) Title: IMPROVEMENTS IN OR RELATING TO SOUND REPRODUCTION SYSTEMS 


ELECTRICAL 
AUDIO IN 



f~\ _ IHOQNCTXUCTKN c,, 



-~<8> 


(57) Abstract 

A stereophonic sound reproduction system is aimed at synthesizing at a multiplicity of points (L) in the listening space an 
auditory effect obtaining at corresponding points (K) in the recording space, to compensate for crosstalk (C 2 i, C 12 ) between the 
loudspeakers, the acoustic response of the listening space, and imperfections in the frequency response of the speaker chan- 
nels Each speaker channel incorporates a digital filter (H) the characteristics of which are ajusted in response to measure- 
ments of the reproduced field. The digital filters are provided by an inverse filter matrix (H) of which the matrix elements 
(H ) are determined by a least squares technique, applied so as to minimise the sum of the squared errors resulting from the 
difference between the actual system impulse response and a desired impulse response. A full bandwidth signal is transmit- 
tprf hv a hvna* s mute for combination with the output signal (9) from the filter, the bypass route including a delay means 
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IMPROVEMENTS IN OR RELATING TO SOUND REPRODUCTION SYSTEMS 
This invention relates to sound reproduction systems. 

5 The invention is particularly, but not exclusively, concerned with 

the stereophonic reproduction of sound whereby signals recorded at a 
multiplicity of points in the recording space such, for example, as 
the notional ear positions of a head, are reproduced in the 
listening space, by being replayed via a multiplicity of speaker 

10 channels, the system being designed with the aim of synthesizing at 
a multiplicity of points in the. listening space an auditory effect 
obtaining at corresponding points in the recording space. 

Imperfections in current two-channel stereophonic reproduction are 
15 seen to arise principally as a result of three causes: 

Firstly, all frequency components of the signal played at each loud- 
speaker channel are reproduced both at the right and at the left 
ears of the listener as though the loudspeakers were the signal 
20 sources of the relevant frequency components, regardless of the 
positions at which these components were created in the recording 
space. 

Secondly, the acoustic response of the listening space provides a 
25 reverberant field which is additive with the reverberant field of 
the recording space. 

Thirdly, "imperfections' 1 in the frequency, response of the speaker 
channels. 


30 


According to cne aspect of the invention, a sound reproduction system 
comprises means for employing a measurement of the reproduced field 
so arranged as to enhance the accuracy of the reproduction system. 


35 
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Preferably the system is arranged to provide stereophonic 
reproduction of sound, whereby signals recorded at a plurality of 
points in the recording space are sought to be reproduced at 
specific points in the listening space by being replayed via a 
5 plurality of speaker channels. 

Each of the channels preferably includes a digital filter the 
characteristics of which are adjusted or set in response to 
measurements of the reproduced field. 

In some situations it may be possible to set once and for all the^ 
filter characteristics but usually the reproduced field will be 
monitored continuously or semi -continuously and the characteristics 
of the digital filter are adjusted in response to an error signal 
from the monitor. 

The filters are preferably adapted to compensate, at least to a 
significant degree, for acoustical cross-talk in the reproduced 
field. 
20 

The filters are preferably provided by an inverse filter matrix of 
which the matrix elements are determined, or have been determined, 
by a least squares technique. 

25 The least squares technique is applied so as to minimise the sum of 
the squared errors resulting from the difference between the actual 
system impulse response and the desired impulse response. 

Any required impulse response may be selected or pre-selected. 
30 Thus, for example, a desired reverberation of the room could be 
chosen. 

In a preferred arrangement the digital filter is preceded by a 
low-pass filter, and a full band-width signal is transmitted by a 
35 bypass route for combination with the output signal from the filter, 
the bypass route including a delay means. 
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The input to the digital filter is preferably down-sampled to reduce 
the sample rate and thereby reduce the amount of processing 
associated with the operation of the filter, the output from the 
filter being interpolated to increase the sample rate. 

5 

According to a second aspect of the invention we provide an 
optimiser unit for interposing between speaker outputs of an 
existing stereophonic reproduction system and the speakers thereof 
to produce a system as claimed in any of the preceding claims, the 
10 optimiser unit having respective optimiser inputs for connection to 

the speaker outputs, respective optimiser outputs for connection to 
the speakers, and an input for receiving a microphone signal from a 
microphone placed in the reproduced field. 

15 One embodiment of the invention comprises a sound reproduction 
system for reproducing from signals x R respectively recorded at K 
points in the recording space, an approximation to these signals at 
L points in the listening space using M loudspeaker channels, where 
K the K recorded signals x R being operated on by an M x K 

20 matrix of compensating inverse digital filters, prior to- their 

transmission via an L x M matrix of transmission paths, the said 
inverse digital filters serving to implement an algorithm operating 
at the inputs to the loudspeaker channels to introduce compensation 
for acoustical cross-talk, acoustical response of the listening 

25 space, and corruptive frequency response characteristics of the 

loudspeaker channels. 

In the limiting case of currently used stereophonic reproduction 
equipment VL = 2, and, in general, L = 2. 

30 


35 


WO 90/00851 


- A- 


!• nACKGR"UND TO THE IKVEXTinN 

The general field of the invention is the stereophonic repro- 
duction of sound in which signals recorded at spaced points in 
the existing sound field such, for example, as the ears of a 
dummy head, ^re replayed via two loudspeakers in a listening 
5 space with the aim of creating a faithful reproduction of the 
original recording space sound at the ears of a listener. 

It has been recognised that imperfections in the reproduction 
arise from" three main causes: 

(1) The signal played via the right channel is reproduced 
10 at both the right and the left ears of the listener. 

Similarly, for the signal played via the left channel. 

(2) The acoustic response of the listening room provides a 
reverberant field which is in addition to the reverb- 
erant field of the space in which the existing re- 

15 cordings were made. 

(3) The frequency response of the loudspeakers used for 
reproduction is imperfect. 
It is the objective of the sound reproduction system described 
here to achieve the "perfect" reproduction of the recorded 
20 signals at the listener's ears. The signals recorded at two 

points in the recording space are to be reproduced essentially at 
exactly two points in the listening space. 

In order to compensate for the three adverse factors mentioned 
above, inverse filters are introduced acting as inputs of the 

25 loudspeaker channels used for reproduction in order to compensate 
both for loudspeaker response and as room response. 
Initial attempts to design such inverse filters for the single 
channel case have been reported by Farnsworth et al (1) and 
Clarkson et al (2). 

30 The formulation proposed includes the effective cancellation of 

undesirable 'acoustical cross talk 1 provided by the unwanted trans- 
mission previously mentioned from the right channel to the left 
ear and vice versa . 

The three problems are thus dealt with in a single formulation 
35 which involves, as made apparent hereinafter, the use of a matrix 
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of digital fillers operating on i he recorded signals prior to 
their transmission via the loudspeaker channels. 
: Vcvi'ii; attempts at addressing certain aspects of the problem 
of stereophonic fidelity have been made using analogue electronics 
5 notably by Atal and Schroeder (o), Damaske (4) and, more recently 
by Sakamoto e_t al^ (5,6). 

In the formulation presented here, however, the use of digital 
techniques is envisaged and the problem can be analysed using a 
least squares approach. 

10 Two novel solutions to the problem are hereinafter presented, 
one involving the design of the inverse filter matrix via a 
,! de terministic" least squares method and the other involving the 
use of a "statistical 11 least squares technique. The second 
approach has the merit that it can be made adaptive and offers 

15 the possibility of computing the filter matrix in situ in a given 
listening space. The proposed system consists of a simple 
addition to existing stereophonic sound reproduction equipment. 
In simple terms the additional filter matrix can be housed in a 
"black box" containing the necessary micro-processors programmed 

20 to implement algorithms hereinafter mentioned so as to operate on 
the two output signals produced by existing equipment. The 
n black box" produces two modified outputs for transmission via 
the loudspeaker system of the existing equipment. 
Additional inputs to the black box are provided by two (or more) 

25 microphones which may be placed at the two (or more) positions 
in the listening space at which the optimal sound reproduction 
is required. These microphones may be removed af'er adaptation 

of the necessary digital filters. 
Further, the approach presented has been generalised such that 

30 it can be applied to any recording technique involving sensing 
the sound field with any number or type of transducers. 
It is thus not only applicable to dummy head or coincident 
microphone techniques but also more recent approaches involving 
the effective detection of the pressure and three components of 

35 particle velocity in the souncJ fieM. 
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In either case the approach presented enables the deduction of the 
inverse filter matrix enabling the production of the best 
approximation (in the least squares sense) to the recorded qualities 
in the listening space. 

5 

The invention will now be further described, by way of example only, 
with reference to the accompanying drawings in which: 

Figure 1 is a diagram showing a stereophonic sound reproduction 
10 system, to illustrate the problem tackled by the invention; 

Figure 2 is a block diagram of a stereophonic sound reproduction 
system prior to the introduction of compensating filters; 

15 Figure 3 is a block diagram of the system of Figure 2 but 

incorporating a matrix H of inverse filters, in accordance with 
the invention, x.^ and x 2 being the recorded signals, and y and 
y^ being the new input signals to the loudspeakers; 

20 Figure 4 is an equivalent block diagram relating to the elements 

of the matrix product H C, and showing the reversal of the order 
of operation of the elements of H and C; 

Figure 5 is a block diagram of a compensated sound reproduction 
25 system in accordance with the invention which incorporates a 

modelling delay of /\ samples to allow for non-minimum phase 
components associated with the transfer functions in the matrix 
C. The desired output signals d^n) and d 2 (n) are now simply 
delayed versions of the recorded signals x^n) and x 2 (n) ; 

30 

Figure 6 is a block diagram of a system in accordance with the 
invention generalised to deal with signals at a plurality K 
points in space and the best least squares reproduction of these 
signals at L points in a listening space, where in general L_>K, 
35 the recorded signals being , the signals output to 

loudpseakers being Y m , and the reproduced signals being "d^ ; 
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Figure 7 is a block diagram equivalent to Figure 4 but 

representing the generation of the signal at the i. 'th point in 
the listening space. 

5 

Figure 8 is a block diagram of a full audio-bandwidth system in 

accordance with the invention in which an equalising filter is 

provided only for low and mid-audio frequencies, but a path is 

provided parallel to that through the equalizing filter to 
IQ provide a full bandwidth path; 

Figure 9 is a more detailed block diagram of a practical high 
sample rate full bandwidth audio system working in accordance 
with the diagram of Figure 8; and 


15 


Figure 10 is a yet more detailed diagram of an optimiser in 
accordance with Figure 9. 


20 


25 


30 
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2. THE INVERSE FILTER MATRIX FOR STEREOPHONIC SOUND REPRODUCTION' 

The basic problem is illustrated in Figure 1. The recorded signals 
x x and x z have been produced at, for example, the ears of a dunray 
head. The objective is to produce precisely the same signals at the 
ears of a listener. These signals are represented by the outputs d x 
5 and d z of microphones placed at, or close to, the ears of the 

listener. Figure 2 shows the transmission system in bloclc diagram form, 
where the transfer functions c% m represent the transmission paths from 
the loudspeakers to the ears of the listener. For perfect sound 
reproduction it is assumed that we require d x = x A and d 2 = x 2 . In . 

10 order to achieve this in principle we must operate on the inputs x A 
and x 2 with a matrix of filters having elements H^jc as shown in 
Figure 3. The presence of the diagonal elements H 21 and B 12 can be 
thought o£ as providing the necessary cancellation of the cross talk _ 
signals present in the listening room as a result of the diagonal 

15 elements C zx and c 12 . in order to find the necessary transfer 
function matrix H we can work in the frequency domain and define the 
vector of signals d as 


:][;:] 


= C y (1) 


20 where the vector of signals y is defined as 

LyJ LH 2i h z JLJ 

Combination of these equations yields 

d = chx (3) 


Thus perfect sound reproduction, i.e., such that d = x, requires 
25 CH = r or H = C"- 1 -. Thus in frequency domain terms the necessary filter 
matrix H is simply the inverse of the matrix C. However, we cannot 
compute an exact inverse of the matrix C due to the presence of non-minimum 
phase components in those transfer functions (it is well known, for 
example, that the room acoustic transfer function contains significant 
30 non-minimum phase elements [7]). it is therefore necessary to take a 
least squares approach to the design of the inverse filter matrix and in 
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order to do this it helps greatly if the block diagram representation of 
the composite system shown in Figure 3 can be rearranged. The necessary 
rearrangement can be explained simply using a little matrix algetora. 
Equation (3) can be expanded in full as 


(4) 


and by expanding the product of the matrix c with the matrix H this 
can be written in the form 


IdJ 


(C iA H iA + C iZ H ZJL ) (C^H^ + C 1Z H ZZ )1 fX^j 

(C^H^ + C ZZ H Z1 ) (c ZJL H iZ + c zz e 2Z )J LxJ 


(5) 


It is also evident that this product can be written as 

x z C i2 | TH, 
x 2 C z J 







10 

4z. 


-XiC zi 



H 


22 


(6) 


where a vector representation of the filter matrix H has now been 
adopted* If we now define the signals 5qcC im as the "filtered 
reference signals" r^m* this equation can be further reduced to 


d,l = 


■121 


■221 


*112 
*212 


'122 


■222 


H 


ii 


H 


21 


4 12 


H 


22 


(7) 


15 The corresponding block diagram of the rearranged system is shown in 
Figure 4, where the input signals x x and - x z are first passed through 
the relevant elements of the transfer function matrix C before being 
passed through the elements of the transfer function matrix H. This 
approach is thus essentially that taken by Elliott et at [7] in deriving 

20 the stochastic gradient algorithm used in problems of active control of 
periodic noise and vibration, and has been extended for use with multiple 
input signals by Nelson et at [8]. This block diagram arrangement can 
now be used to determine the optimai finite impulse response filters 
comprising the matrix H. Two alternative least squares approaches can 
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be taken to determine the filters which give the "best** transmission 
properties of the net system. 
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3, DETERMINISTIC LEAST SQUARES FILTER DESIOf 

in this case we determine the filters comprising H which give the 
best least squares fit to a desired system impulse response. Firstly 
note that this requires an exact knowledge of the transfer functions 
comprising C. Now observe that if the inputs x A and x 2 are both 
5 unit impulses (at time index n = O in discrete time) then the filtered 
reference signals r(n) will be given exactly by the impulse responses 
of the relevant transfer functions in C. Thus the net impulse response 
of the system (the outputs d A and d 2 ) are comprised of convolutions of 
the form 


10 c(j)*h(i) = 


C(O) 

0 

0 

c(l) 

C(0) 

0 

c(2) 

c(l) 

c(0) 

C(3) 

C(2) 

c(l) 

• 

C(3) 

C(2) 

• 

• 

C(3) 


c( J-l ) 

0 C(J-l) 

0 0 C( J-l ) 

o 0 0 

O 0 0 


0 

o 

0 

o 

c(0) 

C(l) 
C(2) 
C(3) 


C(J-1 )J 


rh(o) 

h(i) 
h(2) 
h(3) 


h( i-i )J 


(8) 


where the duration of the impulse responses c(j) is assumed to be J 
samples and the impulse response of the compensating filters h(i) is 
assumed to be of duration I samples. Note it is also implicitly assumed, 
that both c(j) and h(i) are causal impulse responses. It therefore 
15 follows from the rearranged block diagram as expressed in equation (7) 
that the net output sequences d x (n) and d 2 (n) due to unit impulses a 
the inputs at n = 0 can thus be written as a composite sequence vector 
defined by 
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d A <0) 


lc XA < 0 ) o 


IC tx (0) 

I 


<=xx<0) 1C A1 <1) C AA <0) lC iA (l) C AA (0) ic^l) " C i4 <0) 


^ A <*-x)i o c xx (*-i)i o c AA (j-i), o c Ax( j-n 

C iA (0) O JC„(0) O IC Zi (0) O t C xz (0) o 

C^D^C^O) JC„(1) C «(0) |C iX (l> C 4i (0) IC«(1)^ C „(0) 

C*j(l) I . C«<1) j . ^C 1A <1) i . \ Cu (lJ 



(9) 


This equation can be written compactly in matrix form as 


d« = C«h 


(10) 


where the subscript '"a" has been used to denote sequence. We now require 
to be a least squares approximation to a desired impulse d. This 

5 desired response could be, for example, simply a delayed version of the 
input x(n). This will allow the incorporation of an appropriate 
"modelling delay" in the design of the inverse filter and therefore the 
reudction of the minimum mean square error associated with the least 
squares estimate (see Wldrow and Steams [10]). Thus, for example, if a 

10 delay at n = 2 is required, the vector d is given by 


d T = co o i o • . . o : o o l o ... o] 


(ii) 


It should also be noted however that any arbitrary required definition of 
d could be adopted and in principle a desired reverberation of the room 
could be incorporated into the formulation. 

It may also be useful to select the delay in each of the desired signals 
differently for each microphone and this selection could also be 
accomplished adaptively for example to minimise the mean sauare error. 
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We thus seek to minimise the sum of the squared errors resulting from 
the difference between the actual system impulse response and this 
desired impulse response. This cost function can be defined as 

5 ■ I = (d — d s ) T (d - da) (12) 

■ ! 

which, using equation (9), can be written as 

r = (d - Cflh) T (d - Csh) (13) 

This can be expanded to give the quadratic form defined by 

I = d T d - 2d T c fl h + h T Ca T Cflh (14) 

10 This function is minimised by the optimal filter vector of filter 
coefficients defined by 

h 0 = [Cs^sr^d (15) 
with a corresponding minimum mean square error which is given by 


Io = £ T d - d^ho (16) 

15 Thus t ^ e optimal filter vector h Q can be evaluated by direct inversion 
of the matrix [C^Cg]. 

It should be noted 
that exact measurements of the element comprising Cg, i.e., the 
transmission path impulse responses, are required for the computation. 
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4, STATISTICAL LEAST SQUARES FILTER DESIGN 

In this approach the system is assumed to be supplied with a 
spectrally broad "training signal" via the inputs x A and x 2 . The 
resulting outputs can be written (in discrete time) as 

<3i(n) = r^; d 2 (n) = r 2 *h (17) 

5 where the vectors of filtered reference signals are given by 
£x T =tr iJLA T £a2a T £ii2 T 

(18) 

Ez T - £E Zii T r 22i T £ ^Tj 

Each of the component vectors comprising these vectors of filtered 
reference signals are sequences of the form 

10 £*mk T = C**mk< n ) rimkC**- 1 ) **mk<n-2) *imk<n-l+l)] (19) 

The vector h is the composite tap weight vector as defined above in 
- equation (9). The net output can thus be expressed as the composite 
vector 


d(n) = R h 


(20) 


15 where the vector d(n) and the matrix R are defined by 


d(n) = 


d z (n)J 


R = 


Ti 


(21) 


We now seek the optimal filter vector to minimise the time averaged 
error between the actual and desired outputs, in this case the desired 
output is defined as the signal produced by passing the training signal 
20 through a filter having the desired impulse response. Again this could 
simply be a pure delay, A samples, or any other modified form of impulse 
response function. The net block diagram can thus be written in the 
form shown in Figure 5. in this case we seek to minimise the time 
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averaged aura of squared errors defined by 

J = E[(d(n) - d(n)) T (d(n) - d(n))] (22) 

where E denotes the expectation operator and d(n) is the vector 
defined by d T (n) = [d A (n) d 2 (n)]. Thus, J can be written as 

) 

5 J = E[(d(n) - R h) T (d(n) - R n)] (23) 

which reduces to the quadratic form 

J = E[d T (n)d(n)] - 2E[d T (n)R]h + h T E[R T R]n (24) 
This function has the minimum defined by 

Go = CE[R T R]r A E[R T d(n)] . (25) 
10 and the corresponding minimum mean squared error given by 

Jo = ECd T (n)d(n)] - E[d T (n)RJh 0 (26) 


One ^technique is to use Elliott's 
stochastic gradient algorithm where one updates the coefficients of the 
15 composite filter vector h in accordance with 

h(n + 1) = h(n) - aR e(n) (27) 

where a is the convergence coefficient of the algorithm and the 
instantaneous error vector is defined by 

e(n) = d(n) - d(n) < 28 > 

20 The generalisation of the use of this algorithm follows easily from the 
formulation presented here and the description of the algorithm presented 
in reference [12]. The advantage of this technique is that the 
implementation of the algorithm requires only an imperfect knowledge of 
the elements of the matrix c in order to derive the values of the 
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filtered reference signals used in the coefficient update equation given 
by equation (27). m the single frequency case the accuracy with which 
the values of the filters c must be known is well established. More 
vorX, however, is required to establish the necessary requirements for 
the accuracy of the generation of the filtered reference signals in this 


case.— 
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5. EXTENSION TO MULTIPLE CHANNELS 

The above techniques axe easily extended to deal with the case of 
multiple channels where , for example, the signals are recorded at K 
points in the recording space (producing signals x^) and we wish to 
reproduce a test least squares approximation to these signals (d*) at L 
5 points in the listening space using M loudspeaker channels (with outputs 
Ym)* In general it is assumed that L > M > K. Again we operate on the K 
recorded signals with an M x K matrix H of filters prior to their 
transmission via an L x M matrix C of transmission paths. This is 
illustrated in Figure 6. The formulation of this general problem is most 
10 easily accomplished by again using the transfer function reversal 

technique. The approach adopted here is identical to that presented by 
Nelson et al [9]. Thus in matrix form, again working in the frequency 
domain, we write the vector of M loudpseaXer signals and the vector of li 
reproduced signals as 


15 


2 = H x. 


d = C 2 


(29) 


If these are now written in terms of the columns h^ of H and the 
rows c fi of C then we have 


■£L A 


(30) 


The signal produced at the i'th point in the listening space is thus 
20 given by 


&i = = £i T <fijXi + 5z x z UkXk) 


(31) 


Transposing this scalar then yields 


(32) 
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This then demonstrates that the equivalent "reversed transfer function" 
block diagram can be expressed as shown in Figure 7. in discrete time, 
the signal at the I'th point is the listening space can thus be expressed 
in terms of the convolution 


dje(n) = rjTh 


(33) 


where the composite tap weight vector h and the I f th filtered reference 
signal vector rj are defined by 


^ T r T ' T 


T , T T 
E*M± l££iz Elzz • 


T T T T 


(34) 


T T 


T T T T T T T 

• !lMx hzz • • • £tl2 I • * I^jJC hzK • • • ]2hk3 


10 Where each of the component vectors are given by the sequences 


limJc = t^imjcCn) r £nilc (n-l) ... n-l+l)] 


(35) 


The net vector of discrete time signals produced at the 1 points in the 
listening space can now be written as 

15 d(n) = R h (36) 

where the vector d( n ) and the matrix R are now defined as 


d(n) = 




R = 



d*(n) 


m 


• 


a*. 


(37) 


The optimal composite tap weight vector which minimises a cost function 
of the form 


20 


J = E[(d(n) - d(n)) T (d(n) - a(n))] 


(38) 


WO 90/00851 


- 1 9~ 


PCI7GB89/00773 


(where d(n) is now the i^vector of desired impulse responses) is now 
deduced by following exactly the analysis presented in equations (22) to 
(26) above. The adaptive coefficient update equation (27) also remains 
applicable in this general case. — 
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6. PRACTICAL IMPLEMENTATIONS OF THE INVENTIVE SYSTEMS 

Adaptive equalisation and cross-talk suppression of a full 
audio-bandwidth system may be very costly if the equaliser must 

5 operate over the full audio bandwidth. If only low and mid-audio 

frequency equalisation is required then only the frequencies which 
are to be equalised need to be passed through the equalising filter. 
A full bandwidth path must, ' however, be provided for the audio 
signals, the equalisation preferably being achieved by adding the 

10 suitably filtered low frequency components to the total signal so as 
to achieve the desired result at a listening point, as shown in 
Figure 8. 

It should be noted that the required filters for minimum distortion 
15 and cross-talk will result in minimum error and so, in accordance 
with the invention, the error may be used as a criterion for 
adjusting the filters. 

A block diagram of a practical audio response system based on the 
20 system of Figure 8 is shown in Figure 9. 

Full bandwidth audio signals (A) are anti-alias filtered by low pass 
filters (B) and converted into digital form by analogue-to-digital 
converters (C). These steps may be unnecessary, and can therefore be 
25 omitted in a contemplated modification, since digital audio is 
commonly available from compact disc players and digital audio tape 
machines, thus removing the need to digitise analogue audio signals. 

The full audio bandwidth signals (3) are then sent both to the low 
30 frequency optimiser and through a delay (F) to form part of the 
output of the optimiser. The delay (F) is required to ensure that, 
despite any delays that the low frequency components experience 
during their filtering in the optimiser, both high and low frequency 
components may emerge to be added in at the same time. 
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Before being processed by the low frequency optimiser, the full 
audio bandwidth signals (3) are preferably low pass filtered and 
down-sampled so as to reduce the sample rate and the amount of 
processing the optimiser must perform. The low bandwidth signals 

5 (5) are then fed both into the low frequency equaliser (G) which 

consists of an array of adaptive finite impulse response filters, 
and into ^the modelling delay (I). The modelling delay allows the 
filters (G) prior knowledge of the input signals compared with 
the modelling delay outputs (6). This allows the filters (G) to 

10 implement transfer functions with an apparently non-casual part. 

The outputs of the modelling delay are compared with the digital 
representation of the low frequency components of the signals at the 
actual listening points (15). The error (8) between these signals 

15 is then used to drive the adaptive filters (G) so as to minimise 
itself. The outputs of the adaptive filters (9) are then 
interpolated so that the sample rate of these signals may be 
increased to that of the full audio bandwidth signal (11). The 
signals Til) and (10) are then added to produce the output required 

20 for optimum cross-talk suppression and equalisation at the 

listening point (P). 

The digital signals (13) are then converted into analogue form and 
filtered to remove high frequency distortions produced in the 
25 conversion process. The analogue signals are then ready for 
amplification and playing through loudspeakers in the listening 
environment (0). 

It can be seen that this system prefilters and combines the signals 
30 so that cross-talk and amplitude response distortions are removed as 
far as the listener is concerned. 

The block diagram of Figure 9 has been expanded in Figure 10 to give 
a more detailed diagram of a system for stereo reproduction. Most 
35 of this diagram is dealt with in the discussion of the simpler block 
diagram of Figure 9. It is, however, worth discussing the adaptive 
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filters H 1X , h* 21 and H^. 

The parameters of these filters may be found using a modified 
"filtered-x" LMS algorithm. In order to do this a model of the 
5 external electro-acoustic transfer functions C^, C i2' C 21 and C 22 

is required. This may also be found adaptively. 

The filters and H 22 simply equalise the responses of the two 

channels. It should be noted that the equalisation is done in the 
10 time . domain and so the gross phase distortions which accompany 
"graphic", parametric and other analogue filter equalisation- 
techniques will be avoided. 

It should be noted that two sample rates are present in this digital 
15 system. We consider that it may be advantageous to use two 
processors to perform the necessary filtering. The dotted line of 
Figure 9 suggests a preferred divide for the processing - one 
processor handling the optimising while the other takes care of 
delay equalisation, down-sampling, interpolation and up-sampling. 

20 


25 


30 


35 
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CLAIMS 

1. A sound reproduction system characterised by means (H) for 
employing a measurement of the reproduced field so arranged as to 

5 enhance the accuracy of the reproduction system. 

2. A system as claimed in claim 1 characterised in that the system 
is arranged to provide stereophonic reproduction of sound, whereby 
signals recorded at a plurality of points (K) in the recording space 

10 are sought to be reproduced at specific points (L) in the listening 
space by being replayed via a plurality of speaker channels y . ■ 

3. A system as claimed in claim 2 characterised in that each of the 
channels includes a digital filter (H) the characteristics of which 

15 are adjusted or set in response to measurements of the reproduced 
field. 

4. A system as claimed in claim 3 characterised in that the 
reproduced field is monitored continuously or semi-continuous ly and 

20 the characteristics of the digital filter are adjusted in response 

to an error signal from the monitor. 

5. A system as claimed in claim 3 or claim 4 characterised in that 
the filters are adapted to compensate, at least to a significant 

25 degree, for acoustical cross- talk in the- reproduced field, 

6. A system as claimed in claim 5 characterised in that the 
compensation for acoustical cross-talk is continuously or 
semi-continuously reviewed, by monitoring of the reproduced field, 

30 to produce an error signal and any adjustments found desirable are 
automatically made to the relevant filter characteristics in a 
direction to reduce the error signal. 

7. A system as claimed in any one of claims 3 to 6 characterised in 
35 that the filters are adapted to compensate, at least to a 

significant degree, for the acoustical response of the space in 
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which loudspeakers supplied by the respective speaker channels are 
situated. 

8. A system as claimed in any one of claims 3 to 7 characterised in 
5 that the filters are provided by an inverse filter matrix (H) of 

which the matrix elements (H.J are determined, or have been 
determined, by a least squares technique, 

9. A system as claimed in claim 8 characterised in that the least 
10 squares technique is applied so' as to minimise the sum of the 

squared errors resulting from the difference between the actual, 
system impulse response and a desired impulse response. 

10. A system as claimed in claim -9 characterised in that the desired 
15 impulse response is a desired reverberation of a room. 

11. A system as claimed in claim 8 characterised in that the least 
squares technique is applied by considering the system response to a 
spectrally broad input signal, the minimum mean squared error being 

20 determined by use of an equation of the form 

Jo = E [d T (n)d(n) ] - E [ d T (n)R J ho 
the filter coefficients being determined by an equation of the form 

25 

ho = C E C- T -n]" 1 E [£ T ^ (n) 3 

and employing an algorithm of the form 
30 h(n+l) = h(n) - ocRe(n) 

to determine said minimum. 

12. A system as claimed in any one of claims 3 to 11 characterised 
35 in that the digital filter is preceded by a low-pass filter (D), and 

a full band-width signal is transmitted by a bypass route for 
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combination with the output signal from the filter, the bypass route 
including a delay means (F). 

13. A system as claimed in any one of claims 3 to 12 characterised 
5 in that the input to the digital filter is down-sampled (at E) to 

reduce the sample rate and thereby reduce the amount of processing 
associated with the operation of the filter, the output from the 
filter being interpolated (at H, Figure 9) to increase the sample 
rate. 

10 

14. A system as claimed in any one of claims 3 to 13 characterised, 
in that the characteristics of the filter are adjusted utilising a 
modelling delay (at I, Figure 9) which provides a delayed version of 
the input signal to the filter, 

15 

15. An optimiser unit for interposing between speaker outputs of an 
existing stereophonic reproduction system and the speakers thereof 
to produce a system as claimed in any of the preceding claims, the 

20 optimiser unit having respective optimiser inputs for connection to 

the speaker outputs, respective optimiser outputs for connection to 
the speakers, and an input for receiving a microphone signal from a 
microphone placed in the reproduced field. 

25 16. A sound reproduction system as claimed in claim 1, comprising 

inverse digital filter means incorporated in the inputs of several 
loudspeaker channels and operable so as to exhibit a filter 
characteristic conforming to an algorithm such that input signals to 
the said channels are modified in amplitude and phase in a manner 

30 such that at at least one pair of points, spaced apart in the 

listening space by a distance corresponding to the spacing of the 
human ears, the outputs from the loudspeaker channels are mutually 
compensatory, thereby diminishing cross-talk, compensating for the 
reverberant field of the listening space, and compensating also for 

35 the frequency response of the speaker channels. 
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